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Abstract 

i-S^ . Positivity constraints, derived initially assuming parity conservation, for the inclusive reaction of 

the type A(spin 1/2) + _B(spin 1/2) — ^ C + X, where the spins of both initial spin-1/2 particles 
^ |. can be in any possible directions and no polarization is observed in the final state, are generalized 

^ ' to the case of parity violation. By means of a systematic method, we obtain non-trivial bounds 

involving all the spin observables of the reaction and we discuss some relevant physics processes. 
Particularly we discover a non-trivial positivity constraint for the processes, pp — > /Z^ -I- X or 
CN , pp i"^ + X where £^ decayed from /Z'^, which could be checked at the ongoing longitudinal 

' spin program at RHIC. 
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1 Introduction 



Positivity constraints have been widely studied in spin physics to pin down the smallest allowed domain 
^ . for spin observables. This powerful tool has a broad range of applications for exclusive reactions as 
well as for inclusive reactions. It can be used to test the consistency of a given set of available 
measurements and also the validity of specific dynamical assumptions in theoretical models. Different 
methods can be used to establish these constraints and many interesting cases have been presented in 
a recent review article [T]. 

In the present paper, we will focus on the single particle inclusive production in polarized hadronic 
collisions, ^(spin l/2)-|-i?(spin 1/2) — )■ C+X, where only initial spins are observed and no polarization 
for the final state particle C is measured. If parity is conserved, this reaction is fully described in 
terms of eight independent spin observables. The positivity constraints for these parity- conserving 
observables have been derived in Ref. [2J. However, if parity is not conserved, there are twice as 
many spin observables and the positivity constraints become much more involved. Nevertheless, such 
a case has been partially studied in Ref. p]. The derived positivity constraint has been further used 
to put non-trivial bounds on several Sivers functions entering the theoretical description of single 
transverse spin asymmetry for various processes [IHS]. However, the increasing complexity for parity- 
violation case requires a systematic method to obtain all the positivity constraints. This will be the 
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aim of this paper, i.e., deriving all the positivity constraints for the most general cases - including 
both parity- conserving and parity-violating spin observables. 

The longitudinal program at Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory is currently under successful running, aiming to pin down the polarized antiquark dis- 
tribution in the proton |10[|llj . The reason that boson production could provide unique and 
clean access to the individual antiquark polarizations is due to the maximal violation of parity in the 
elementary Wqq' vertex [12^il3j. Because of the parity-violation nature of this process, pp — )• + X 
ov pp ^ £^ -\- X where decayed from W^, our newly derived positivity constraints could have 
interesting and nontrivial implications for the spin asymmetries measured in these processes. Spin 
asymmetries for the process pp ^ + X can be also measured, although the production rate is lower 
than for production. 

The remainder of this paper is organized as follows: in the next section we review the derivation 
of positivity constraints for the parity-conservation case, to introduce the notation and also to update 
the method which could be easily used for the parity-violation case. In section 3, we will derive 
all the general positivity constraints, which involve both parity-conserving and parity-violating spin 
observables, and are classified according to different degrees, linear, quadratic, cubic and quartic. 
We then give one phenomenological example of our derived positivity bound, when applying to the 
inclusive W^/Z^ production in longitudinal pp scatterings. We summarize our results in section 5. 

2 Positivity constraints for parity-conserving processes 

For single particle inclusive production: ^(spin 1/2) + i?(spin 1/2) C + X with spin vectors Pa 
and Pb for initial particles A and B, respectively, the spin-dependent cross section a{Pa,Pb) can be 
defined through the cross section matrix M and the spin density matrix p 

a{Pa,Pb) = Ti{Mp), (1) 

where p = pa ® Pb '^^ the spin density matrix with pi = {I2 + Pi ■ i?)/2, i = a,b. Here a = {ax, o'y, az) 
stands for the three 2x2 Pauli matrices and I2 is the 2x2 unit matrix, thus p being the direct product 
of Pa and Pb will be 4 x 4 matrix. We will now study the parametrization of the 4x4 cross section 
matrix M. We first review the existing study for the parity-conserving case. The purpose is twofold: 
first to set up notation; and second, to update the method for deriving the positivity constraints, 
which could also be used in the parity-violation case. 

For the parity- conserving process, M could be parametrized in the following way 

M = (To [h + AaNCTay <^ h + Ab^h CFby + ANNC^ay ® CTby + A^aaz tX" (Jbz 
+ AsS(Tax (8) CTbx + ALSCTaz CTbx + AslCT 

Here I4 is the 4x4 unit matrix and ao stands for the spin-averaged cross section. In other words, 
for a parity-conserving process, there are eight independent spin-dependent observables [Il[2l[T3]: the 
unpolarized cross section ao, two single transverse spin asymmetries AaN and AbN, and five double 
spin asymmetries Ann, All, Ass, Als, and Asl- Here the subscript L, N, S represents the unit 
vectors along the spin directions of initial particles A and B. Specifically in the center-of-mass system 
of A and B, L, N, S are along the incident momentum, along the normal to the scattering plane 
which contains A, B and C, and along N x L, respectively. The expression in Eq. ([2]) is fully justified, 
since we have explicitly 

a{Pa,Pb) = ao[l + AaNPay + AbNPby+ ANNPayPby + ALLPazPbz 

+AssPaxPbx + ALsPazPbx + AsLPaxPbz] ■ (3) 
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The positivity constraints for the parity-conserving process have been derived in [2]. The crucial 
point is that the cross section matrix M is a Hermitian and positive matrix. The necessary and 
sufficient condition for a Hermitian matrix to be positive is that all its eigenvalues are positive. It 
is important to emphasize here that the eigenvalues of a matrix are independent of the basis where 
it is written. In other words, no matter what basis one chooses to express the cross section matrix 
M, one should obtain the same set of eigenvalues, from which one obtains the same set of positivity 
bounds. Since this is the necessary and sufficient condition, the bounds derived from this should be 
the strongest constraints for the spin observables. Using this fact about the eigenvalues, we thus could 
choose any convenient basis such that the matrix has a simple form and thus the eigenvalues could 
be easily derived from there. For example, the original positivity bounds are derived by choosing the 
transverse basis where ay is diagonal, in which the matrix elements Mij of the cross section matrix M 
are given by 

Mil = (1 + Ann) + (AaN + Am) , 
M22 = (1 - Ann) + {AaN - Aiv) , 
M33 = (1 - Ann) - {AaN - Aiv) , 

M44 = (1 + Ann) - {AaN + Am) , (4) 

Mi4 = M|i = All - Ass - i {Asl + Als) , 

M23 = M3*2 = All + Ass - i {Asl - Als) , 

M12 = M21 = Mi3 = M31 = M24 = M42 = M34 = M43 = 0, 

i.e., half of the matrix elements vanish in this basis. The eigenvalues Aj can be easily obtained and 
are given by 

Ai,2 = 1 + Ann ± V{AaN + Am)^ + {All - Ass)^ + {Als + Asl)^ (5) 
A3,4 = 1-Ann± V{AaN - Am)^ + {All + Ass)^ + {Als - Asl)^ (6) 

Then from all the eigenvalues A, > 0, we have the following strongest positivity constraints 

(1 ± Ann)^ > {AaN ± Am)'' + {All T Ass)"" + {Als ± ^sl)', (7) 

which are exactly the same as those derived in 

It is nice to be able to derive the eigenvalues of the cross section matrix and thus to obtain the 
strongest positivity bounds directly. However, sometimes the eigenvalues turn out to be very difficult 
to find, which is exactly the situation we face for the parity-violation case in the next section. In 
such situations, one could derive a complete set of necessary constraints which as a whole forms the 
sufficient condition for positivity. This is the strategy of J.J. Sylvester, so-called Sylvester's criterion, 
which states that the necessary and sufficient condition for a Hermitian matrix to be semi-positive is 
that all its principal minors have to be non-negative. This is equivalent to say for the 4x4 matrix M 
that: 

• All the diagonal matrix elements Ma > 0; 

• The elements satisfy: MaMjj > |Mjjp; 

• The determinant of any 3x3 matrix formed by removing from M its i^^ row and i*^ column are 
non-negative, where i = 1,2,3, or 4; 

• The determinant of the matrix M itself is non-negative. 



3 



These conditions will enable us to derive all the general positivity constraints, which are of different 
degree, linear, quadratic, cubic and quartic, respectively, and as a whole they form the necessary and 
sufficient condition for the positivity of M. 

Let's use Sylvester's criterion to revisit the positivity bounds for parity-conserving case. First let's 
work again in the transverse basis. With the expression of M in Eq. ^ in hand, from Ma > we 
immediately derive 

l±ANN>\AaN±AbN\. (8) 

On the other hand, from MaMjj > |Mjjp, we further derive 

(1 ± Ann? > (AaN ± AtN? + {All T Ass? + {Asl ± Als?, (9) 

It is also easy to find that going to the even higher order principal minors does not give any further 
constraints. Thus we have derived the equivalent results from a slightly different method. 

We would also like to make some comments on the positivity bound in Eq. ([8]). Although this 
is a weaker bound and could be deduced from Eq. ([7]), its derivation is much simpler in the sense 
that one only deals with the diagonal matrix elements without looking for the eigenvalues which need 
extra (and difficult) mathematical work. On the other hand, these kinds of linear bounds are also 
very useful in the phenomenological studies because of the limited accessibility and accuracy of the 
spin asymmetries in the experimental measurements. 

It is also important to realize that Sylvester's criterion holds true in any basis one chooses to 
express the cross section matrix M. Thus choosing a different basis, one could derive a different 
set of positivity bounds. Even though these different sets are equivalent to each other according to 
Sylvester's criterion, one might obtain some interesting positivity bounds by choosing a convenient 
basis, which could be difficult to find directly in another basis. 

For example, if we choose the helicity basis where cr^ is diagonal, the cross section matrix M 
becomes 



Mil 


= M44 


= 1 + All 


M22 


= M33 


= 1-All 


Mi4 


= M4i 


= Ass — Ann 


M23 


= M32 


= Ass + Ann 


M12 


= M*, 


= Als - iAbN 


Mi3 


= ^3*1 


= AsL + iAaN 


M24 


= M|2 


= -AsL - iAaN 


M34 


= M|3 


= -Als - iAbN 



(10) 



Now Ma > leads to the trivial bounds 1 it All > 0. However, using MuMjj > |Mijp leads to some 
new positivity bounds: 

1-AIl > AIs + AIn, (11) 
1-AIl > A^^l + AIn, (12) 
{l±ALLf > {Ass T Ann f, (13) 

which are not very easy to derive in the transverse basis. These are actually weaker bounds compared 
with the strongest bounds in Eq. ([7]). It requires some extra work to derive them in the transverse 
basis: if one uses the following formula in Eq. ([7]) 

^Ja'l+al~\ h > |ai + 02 H h , (14) 
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one could derive all the inequalities (jlip - ()13p . One could continue to study the higher order principal 
minors and thus to obtain the cubic and quartic constraints. Then together with the linear and 
quadratic ones, they form the complete set of positivity bounds which should be equivalent to the 
bounds in Eq. ([7]) according to Sylvester's criterion. 

3 All general positivity constraints 

We now study the general positivity constraints, which involve also the parity-violating asymmetries 
where one will have sixteen independent spin-dependent observables [3j. Besides those in the parity- 
conserving processes, one has four additional single spin asymmetries Aai, -A-bL, -^aS and A^s, and 
four additional double spin asymmetries ^lat, ^atl, A^s and Asn- In this case, the more general 
cross section matrix M can be parametrized as 

M = O-Q [h + AaNCTay <^ h + Amh CTby + AN^CTay <^ CTby + AiLf^az <X) (Tbz 
+AsSCrax O (Jbx + ALSCTaz O CTbx + AsLCT ax O'bz] 
+C70 [AaLCTaz h + Ablh ® CFbz + AaSCTax ®h + Absh ® CFbx 

+ALNO-az (X" CTby + ANL(^ay (X" (Tbz + ANSO^ay (X) (Tbx + AsN(Tax (^by] , (15) 

which is fully justified since one has 

(T{Pa,Pb) = CTqII + AaNPay + AbNPb-y + A^NPayPby + ALLPazPbz 
+AssPaxPax + AisPazPbx + AsiPaxPbz] 
+aQ[AaLPaz + AblPbz + AasPax + ^fesAx 

+ALNPazPby + ANlPayPbz + ANsPayPbx + AgNPaxPby] ■ (16) 

Again to study the positivity constraints, we need to express the cross section matrix M in a basis and 
then study the necessary and sufficient conditions for M to be positive. Because of the complexity of 
the matrix M, it becomes difficult to obtain the eigenvalues. Instead we will use Sylvester's criterion. 

First, we will still work in the transverse basis where ay is diagonal, in which the matrix elements 
of M are given by 

Mil = (1 + Ann) + {AaN + Atv) 
M22 = (1 - Ann) + {AaN - Atv) 
M33 = (1 - Ann) - {AaN - Atv) 
M44 = (1 + Ann) - {AaN + Atv) 
Mi4 = = All -Ass-i {Asl + Als) 

M23 = M^^ = ALL + Ass-i{AsL- Als) ^ ' 

M12 = M*^ = AbL + Anl - i {Abs + Ans) 

Mi3 = Mg*! = AaL + Aln - i {AaS + Asn) 

M24 = = AaL - Aln - i {Aas - Asn) 
M34 = MX3 = AbL - Anl - i {Abs - Ans) 

Start with the linear positivity bounds. From Ma > 0, we immediately obtain 

1±Ann >\AaN±AbN\- (18) 

The quadratic bounds are derived from MaMjj > |Mjjp. Especially for = {1,4}, and {2,3}, 

one obtains 

(1 ± Ann)^ > {AaN ± AbN)^ + {All T Ass)^ + {Asl ± ^ls)'- (19) 
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Both Eqs. (jlSp and ()19p are exactly the same as those in Eqs. ^ and ([9]), which are derived for the 
parity-conserving case. Since all the asymmetries involved in these inequalities are parity-conserving 
ones, it is not surprising that these positivity bounds are preserved. For the case when {i,j} = 
{1, 2}, {1,3}, {2, 4}, and {3, 4}, we obtain four extra positivity bounds 

il± AaN? > {AbN±ANNf + {AbL±ANLf + iAbS±ANsf, (20) 
{l±AhNf > {AaN±ANNf + {AaL±ALNf + {AaS±AsNf, (21) 

which involve both parity-conserving and parity-violating asymmetries. 

The cubic bounds are derived from the determinants of all the 3x3 matrix formed by removing 
from M its i^^ row and column, which has to be non-negative. If we remove the 4*^^ row and column, 
i.e., keep the row 1, 2, 3 and column 1, 2, 3, this newly formed matrix denoted as M123 is given by 

/ AaN + AbN + Ann + 1 + Anl - i{Abs + Ans) AaL + Aln - i{Aas + Asn) \ 

Mi23= At,L+ANL+i{Abs + ANs) AaN- AbN- Ann + I All + Ass - iiAgL - Als) (22) 

\ AaL + ALN + iiAas + AsN) All + Ass + iiAsL - Als) -AaN + AbN - Ann + ^ J 

Then the determinant IM123I > leads to the following inequality 

IM123I = (1 + Ann + AaN + AbN)[{l - Ann)^ - {AaN - AbNf - {All + Assf - {Asl - Als?] 
-{AaN - AbN - Ann + 1)[(Al + Anl? + (As + Ans?] 
-{-AaN + AbN - Ann + 1)[(Al + Aln? + {Aas + Asn?] 
+2{All + Ass)[{AbL + ANL){AaL + Aln) + (As + ^ivs)(As + Asn)] 
+2{Als - AsL)[{Abs + ANs){AaL + Aln) - (Al + ^iVL)(As + Asn)] > (23) 

Likewise, from |-/\fi24|, I-M1341, 1^^234 1 ^ 0, we have 

IM124I = (1 - Ann + AaN - Ajv)[(1 + Ann)^ - {AaN + Aiv)^ - {All - Ass)^ - {Asl + Als)^] 
-{1 + Ann - AaN - A7v)[(Al + Anl)"^ + (As + ^tvs)^] 
-(1 + Ann + AaN + Aiv)[(AL - Aln)^ + (As - Asn)^] 
+2(Al - ALN)[{AbL + Anl){All - Ass) + (As + Ans){Asl + As)] 
+2(As - AsN)[{AbL + Anl){Asl + Als) - (As + Ans){All - As)] > (24) 

IM134I = (1 - Ann - AaN + Ajv)[(1 + ^^vtv)^ - (Aiv + Aiv)^ - (Al - As)^ - (Al + As)^] 
-(1 + Ann - AaN - Aiv)[(AL + Aiv)^ + (As + Aiv)^] 
-(1 + Ann + AaN + A7v)[(Al - Anl)^ + (As - Ans)^] 

+2{AbL - ANL)[{AaL + Aln){AlL - As) + (As + A7v)(Al + As)] 

+2(As - ANs)[{AaL + Aln){Asl + Als) - (As + A7v)(Al - As)] > (25) 

IM234I = (1 + Ann - AaN - Aiv)[(l - Ann)^ - {AaN - Aiv)^ - {All + As)^ - {Asl - Als)^] 
-{1 - Ann - AaN + Aiv)[(AL - Aiv)^ + (As - Aiv)^] 
-(1 - Ann + AaN - A7v)[(Al + Anl)^ + (As + Ans)^] 
+2{AbL - Anl)[{All + As)(Al - Atv) + (Al - As) (As - Atv)] 
+2(As - Ans)[{All + As) (As - Asn) - {Asl - As)(Al - Aln)] > (26) 

Finally the quartic bounds are given by the determinant of M itself \M\ > 0. With the matrix 
elements given in Eq. (jl7p . it is easy to write down the determinant of M, which we do not write out 
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explicitly here. The linear, quadratic, cubic and quartic bounds form the complete set of positivity 
bounds according to Sylvester's criterion. We have found that all the bounds derived from the parity- 
conserving case are preserved in this more general case. 

As we have emphasized in the previous section, we could also derive a different set of positivity 
bounds by choosing a different basis. These new bounds might also be useful and interesting for the 
phenomenological studies. Let us study the positivity bounds by choosing the helicity basis, in which 
the explicit form of M is given by 

Mil 

M22 
M33 
M44 

Mi4 

M23 
M12 

Mi3 

M24 
M34 

Now from Ma > 0, we have 



This is a very interesting positivity bound. Particularly it involves both parity-conserving {All) and 
parity- violating (Al) asymmetries. Since this bound involves only the longitudinal asymmetries, it 
might be very useful and relevant to the ongoing longitudinal W^/Z^ program at RHIC. 
Then from MuMjj > |Afjjp, one could derive the following quadratic bounds 

{l±ALLf > {AaL±AhLf + iAssTANNf + {ANS±AsNf, (29) 
{l±AaL? > {AbL±ALLf + iAbs±ALS? + iAbN±ALNf, (30) 
{l±AbLf > {AaL±ALLf + {AaS±AsL? + {AaN±ANLf. (31) 

These bounds are stronger than the linear ones. For example, specifically from the bounds in Eq. (j29p . 
one could deduce those in Eq. (j28p . One could continue to study the bounds in even higher order, and 
the procedure is straightforward. We decide to stop here, and turn to discuss some phenomenological 
applications. 

4 Phenomenological example: W^/Z^ production at RHIC 

In the last two sections, we have derived quite a few positivity constraints which involve both parity- 
conserving and parity- violating spin asymmetries. They could have broad applications in testing the 
consistency of the experimental measurements, or studying the validity of the theoretical models. 
However, due to the experimental limited accessibility and accuracy, only a few could be reachable in 
the near future. In this section, we study one such example: the positivity bound (j28p in the W^/Z^ 
or £^ production in pp collisions, pp — ?• W^/Z^ + X or pp ^ £^ + X where £^ decayed from W^/Z^. 
It is important to realize that for identical initial particles scattering, one has 

AaLiy) = AbLi-y), (32) 



= (1 + All) + (AaL + Al) 
= (1 - All) + {AaL - A^l) 
= (1 - All) - {AaL - A^l) 
= (1 + All) - {AaL + Al) 





= Ass 


— Ann - 


- i {Ans + AsN 


= M3*2 


= Ass 


+ Ann - 


- i {Ans - Asn 


= M|i 


= Abs 


+ Als - 


i {AbN + Aln) 


= M*, 


= AaS 


+ AsL - 


i {AaN + Anl) 


= Ml^ 


= AaS 


-AsL - 


i {AaN - Anl) 


= Mi, 


= Abs 


-Als - 


i {AbN - Aln) 



(27) 



l±ALL>\AaL±AbL\. (28) 



7 



where y is the rapidity of the final-state particle. Thus Eq. ()28p becomes 



l±ALL{y)>\AL{y)±AL{-y)\. (33) 

These bounds have a very simple form, and should be very interesting to test them in RHIC experi- 
ments. 

To have an idea, let us check whether these bounds are satisfied for or production in 
longitudinal pp collisions, pp ^ /Z^ + X. In perturbative QCD formalism, at leading-order and 
restricting to only up and down quarks, one has the following simple expressions for the single spin 
asymmetries [ElIiailSHIZ! 

A^^^iy) = + ^d{xa)u{xb) ^^^^ 

u{Xa)d{xh) + d{Xa)u{xh) 
jiW- ,y\ = -^d{Xa)u{xi,) + /\u{Xa)d{Xb) 
^ d{Xa)u{Xb) +u{Xa)d{xh) 

zo _ T.gi'^Vqaq) [-Aq{xa)q{xb) + Aq{xa)q{xb)] 

^ " Eqivl + ^l)[QMqixb) + q{xa)qixb)] ' ^ ' 

and for the double spin asymmetries 

A^~^{y) - Au{xa)Ad{xb) + Ad{xa)Au{xb) ^^^^ 

u{Xa)d{Xb) + d{Xa)u{xb) 

.W-( \ _ Ad{xa)Au{xb) + Au{xa)Ad{xb) , . 

" d{Xa)u{Xb)+uiXa)d{Xb) ' ^ ^ 

^ EgK' + Qg) [Aq{xa)Aq{xb) + Aqixa)Aq{xb)] 
"-"-^^^ " Eqiv'q+a^q)[qixamxb) + qixa)q{xb)] ' ^ ' 

where Aq{x) and q{x) are the helicity distribution and unpolarized parton distribution function, 
respectively. Vg and are the vector and axial couplings of the Z^ boson to the quark. Xa^b are the 
parton momentum fractions given by 

Xa = mq/ \/s e^, Xb = mq/ ^fs e"^, (40) 

with mq , y the mass and rapidity of the W (or Z) boson and ^/s the center-of-mass energy. 

To estimate these asymmetries numerically, we choose BBS2001 polarized and unpolarized parton 
distribution functions based on a statistical approach |18j . At ^/s = 500 GeV, our calculations plotted 
as a function of rapidity y are shown in Fig. [T]for (left) and W~ (right), and in Fig. [2]for Z^ boson, 
respectively. The solid curves are the single longitudinal spin asymmetry Al, the dashed curves are 
the double longitudinal spin asymmetry All, while the dotted curves are the following combination 

l + ALL{y)-\AL{y) + AL{-y)l (41) 



which must be positive according to the positivity bounds in Eq. (I33p . As we can see from both plots 
that even though \AL{y) + AL{—y)\ could become quite sizable, it is still smaller than 1 + ALLiy), 
thus the bound is satisfied in this leading order calculation for both and Z^. The other bound 
1 — ALLiy) ~ \ALiy) — AL{—y)\ could become even more sizable and is also satisfied within BBS2001 
parameterization. 

We have also checked other popular parametrizations of polarized parton distribution functions, 
to see whether they satisfy our bounds. Particularly we have checked GRSV2000 [19j, AAC2008 [20j, 
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Figure 1: Longitudinal asymmetries are plotted as a function of rapidity y of the W boson in collisions: 
(left) and W~ (right). The solid curves are the single longitudinal spin asymmetry Al, the dashed curves 
are the double longitudinal spin asymmetry All, and the dotted curves are the combination of 1 + All{]j) — 
\AL{y) + Al{-v)\. 




-1.5 -1 -0.5 0.5 1 1.5 

y 

Figure 2: Same as Fig. [U but for Z° boson. 



DSSV [2Tl[22], and LSS2010 [23]. For the unpolarized parton distribution functions, we use exactly 
the same set as the one when the global fitting was performed. That is, we use GRV98 [23] for both 
GRSV2000 and AAC2008, while MRST2002 [25] for both DSSV and LSS2010. It turns out that both 
GRSV2000 and AAC2008 satisfy our bounds, while both DSSV and LSS2010 could have violation at 
large rapidity |y|. The violations are shown in Fig.O in which we plot l+Aii{y) — \AL{y)+AL{—y)\ as 
a function of rapidity y for VF"*", W~ , and Z^. Since for identical incoming hadrons ALL{—y) = Aii{y), 
the combination 1 + Aii{y) — \AL{y) + AL{—y)\ is symmetric under y o —y and we thus only plot 
for positive y. We immediately find that when rapidity becomes large y > 1.5 where Xa ^ 0.7 and 
Xh ^ 0.04 (or vice verse for y < —1.5), the combination 1 + Aii{y) — \AL{y) + AL{—y)\ could become 
negative for both DSSV and LSS2010, and for all W^/Z^ bosons. Since our bounds are very general, 
coming from the positivity conditions of the cross section matrix, they should always be satisfied. 
Thus our newly derived bounds, though very simple, immediately put nontrivial constraints on the 
parametrizations of the polarized parton distributions, for both DSSV and LSS2010. 

According to next-to-leading order calculations in |16[ll7j. the QCD radiative corrections for the 
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0.5 1 1.5 0.5 1 1.5 

y y 

Figure 3: Asymmetry 1 + AiLiu) — |^L(y) + y)| are plotted as a function of rapidity y for the 
parametrization of DSSV (left) and LSS2010 [23] (right). The sohd, dashed and dotted lines 

are for W~^, W^, and Z^, respectively. 

asymmetries are small. We thus expect our findings and conclusions for the specific parametrizations 
of polarized parton distribution functions will not alter in higher-order QCD calculations, for both 
W^/Z^, and for the lept ons decayed from them. It will be very interesting to check our bounds in 
the forward or backward regions (when \y\ becomes large) in the experiments. 

We notice that there are already published data for single spin asymmetry Al at mid-rapidity 
y = from both STAR [26] and PHENIX [27], for the leptons decayed from W'^/Z'^ bosons. 

STAR: A'l = -0.27 ± 0.10, Af = 0.14 ± 0.19 (42) 
PHENIX: = -OMtofi, Al = 0.88+°;^^ (43) 

On the other hand, at y = our bounds becomes 

1 + All{0)>2\Al{0)\. (44) 

At the same time All(0) is very small (a few percents) for the current parametrization of polarized 
parton distribution functions. If one takes ^^^(0) ~ in Eq. (I44p . we obtain 

I^l(0)| < ^. (45) 

Comparing with both STAR and PHENIX data, we immediately find that STAR data is consistent 
with our bounds, while PHENIX data (central value) is certainly out of the bound. Of course, so far 
the data has very large uncertainty at the moment. Nonetheless, it will be important to check all our 
bounds in the future experiments. We look forward to the future experimental measurements to test 
these bounds. 

5 Summary 

We derive all the general positivity constraints for the spin observables in the single particle inclusive 
production, A(spin 1/2) -|- i?(spin 1/2) C + X , where the spins of both initial spin-1/2 particles 
can be in any possible directions and no polarization is observed in the final state. By means of a 
systematic method, we generalize the previous positivity constraints derived for the parity-conserving 
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processes to the most general processes including also the parity- violating ones. We find that the 
positivity constraints involving only the parity-conserving asymmetries are preserved in the parity- 
violating case. 

With the help of Sylvester's criterion, we derive all the general positivity constraints, which are of 
different degree, linear, quadratic, cubic and quartic, respectively. These constraints form a complete 
set of necessary and sufficient condition for positivity. As a special example, we discover a very 
interesting non-trivial bound for the parity-conserving and parity-violating longitudinal asymmetries 
All and Al- This could be relevant to the processes, pp — W^/Z^ + X or pp ^ £^ + X where £^ 
decayed from W^/Z^ in longitudinal scatterings, which are currently under active investigation at 
RHIC. 

The positivity constraints derived here could have broad applications in testing the consistency of 
the experimental measurements, or studying the validity of the theoretical models. We look forward 
to the future experimental data to test these positivity bounds. 

Before closing, let us mention another possible application of our results in polarized deep inelastic 
scattering at high energies where weak interactions contributions, both neutral and charged current 
processes are taken into account, as well as the parity violating polarized nucleon structure functions. 
The explicit expressions of the charged lepton asymmetry AeL{y) and the proton asymmetry ApL{y) 
can be obtained from Ref. [28] and our bounds can be used to put new constraints on the helicity 
distributions. These parity-violating asymmetries are also expected to be modified by effects beyond 
the Standard Model, e.g. scalar and vector leptoquarks [29], and our bounds could be used to put 
further constraints on existing models. 

Acknowledgments: We thank Andreas Metz and Werner Vogelsang for discussions. We are grateful 
to RIKEN, Brookhaven National Laboratory, and the U.S. Department of Energy (Contract No. DE- 
AC02-98CH10886) for supporting this work. 

References 

[1] X. Artru, M. Elchikh, J. M. Richard, J. Soffer and O. V. Teryaev, Phys. Rept. 470, 1 (2009) 
|arXiv:n8n2.0lMI [hep-ph]]. 

[2] J. Soffer, Phys. Rev. Lett. 91, 092005 (2003) [ arXiv:hep-ph/0305222] . 

[3] Z. B. Kang and J. Soffer, Phys. Lett. B 695, 275 (2011) [a rXiv: 1003.4913] [hep-ph]]. 

[4] Z. B. Kang and J. W. Qiu, Phys. Rev. Lett. 103, 172001 (2009) [arXiv: 090373629] [hep-ph]]. 

[5] Z. B. Kang and J. W. Qiu, Phys. Rev. D 81, 054020 (2010) [arXiv:0912.1319l [hep-ph]]. 

[6] M. Anselmino et ai, Eur. Phys. J. A 39, 89 (2009) [arXiv: 0805. 2677 [hep-ph]]. 

[7] M. Anselmino, M. Boglione, U. D'Alesio, S. Mefis, F. Murgia and A. Prokudin, Phys. Rev. D 79, 
054010 (2009) [arXiv:0901.3078 [hep-ph]]. 

[8] L. Gamberg and Z. B. Kang, Phys. Lett. B 696, 109 (2011) [arXiv:1009.193 6' [hep-ph]]. 

[9] Z. B. Kang, J. W. Qiu, W. Vogelsang and F. Yuan. [a?Xw:1103.1591l [hep-ph]. 

[10] G. Bunce, N. Saito, J Soffer and W. Vogelsang, Ann. Rev. Nucl. Part. Sci. 50, 525 (2000) 
[arXiv:hep-ph/ 00072]¥| . 



11 



[11] G. Bunce et al, "Plans for the RHIC Spin Physics Program", June 2008, see: 
spin . riken . bnl . gov/ rsc. 

[12] C. Bourrely and J. Soffer, Phys. Lett. B 314, 132 (1993). 

[13] C. Bourrely and J. Soffer, Nucl. Phys. B 423, 329 (1994) [arXiv:h ep-ph/9405250| . 

[14] G. R. Goldstein and J. F. Owens, Nucl. Phys. B 103, 145 (1976). 

[15] T. Gehrmann, Nucl. Phys. B 534, 21 (1998) [arXi v:hep-ph/9710508] . 

[16] D. de Florian and W. Vogelsang, Phys. Rev. D 81, 094020 (2010) [ arXiv: 1003.4533] [hep-ph]]. 

[17] C. von Arx and T. Gehrmann, lax Xiv: 1103. 1465] [hep-ph]. 

[18] C. Bourrely, J. Soffer and F. Buccella, Eur. Phys. J. C 23, 487 (2002) [ arXiv:hep-ph/0109160] . 

[19] M. Gluck, E. Reya, M. Stratmann and W. Vogelsang, Phys. Rev. D 63, 094005 (2001) 
[arXiv:hep-ph/0011215 |. 

[20] M. Hirai and S. Kumano [Asymmetry Analysis Cohaboration], Nucl. Phys. B 813, 106 (2009) 
^arXiv:0808.04 13 [hep-ph]]. 

[21] D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. Lett. 101, 072001 (2008) 
^arXiv:0804.0422 [hep-ph]]. 

[22] D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. D 80, 034030 (2009) 
[arXiv:0904.38 21 [hep-ph]]. 

[23] E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 82, 114018 (2010) [arXiv: 1010.0574] 
[hep-ph]]. 

[24] M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 (1998) |arXiv:hep-ph/9806404| . 

[25] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28, 455 (2003) 
[arXiv:hep-ph/0211080] . 

[26] M. M. Aggarwal et al. [STAR Collaboration], Phys. Rev. Lett. 106, 062002 (2011) 
|arXiv:1009.0326i [hep-ex]]. 

[27] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 106, 062001 (2011) [arXiv: 1009.0505] 
[hep-ex]]. 

[28] M. Anselmino, P. Gambino and J. Kalinowski, Z. Phys. C 64, 267 (1994) |arXiv:hep-ph/9401264] . 
[29] P. Taxil, E. Tugcu and J. M. Virey, Eur. Phys. J. C 14, 165 (2000) |arXiv:hep-ph/9912272] . 



12 



